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We analyze the achievable precision for single-qubit gates that are based on Raman transitions between two 
near-degenerate ground states via a virtually excited state. In particular, we study the errors due to non-perfect 
adiabaticity and due to spontaneous emission from the excited state. For the case of non-adabaticity, we calculate 
the error as a function of the dimensionless parameter \ = At, where A is the detuning of the Raman beams 
and r is the gate time. For the case of spontaneous emission, we give an analytical argument that the gate errors 
are approximately equal to A7/A, where A is the rotation angle of the one-qubit gate and 7 is the spontaneous 
decay rate, and we show numerically that this estimate holds to good approximation. 
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I. INTRODUCTION 

The ability to perform arbitrary unitary transformations on 
individual qubits is very important in the context of quantum 
computation 1 1]. From the point of view of decoherence, it 
is often advantageous to use near-degenerate ground states of 
a given system as qubits, e.g. different hyperfine levels in a 
trapped ion or atom, or the spin states of an individual ex- 
cess electron in a quantum dot. In such situations, the use of 
optical Raman transitions can be an attractive approach for re- 
alizing single-qubit operations, cf. Ref. |2] for experiments 
with trapped ions, Ref. j3| for a proposal involving atoms in 
an optical lattice, and Ref. 0] for a proposal with spins in 
quantum dots. 

A concrete procedure for realizing arbitrary single-spin op- 
erations via Raman transitions has recently been proposed in 
the context of single spins in quantum dots @|- Variations 
involving different excited states of the quantum dots (light- 
hole excitons instead of heavy-hole excitons) were discussed 
in 0] and 01 . Let us note that there are several proposals on 
how to realize optically controlled two-qubitgates between in- 
dividual spins in different quantum dots |J,|a,l3.lMl- Schemes 
for qubit measurement have also been proposed IsSlMI- 

In Ref. 0] it was suggested to perform the Raman operation 
in an adiabatic fashion, in order to minimize the population in 
the excited state, and thus gate errors due to the decoherence 
of that state, which is in general much faster than the deco- 
herence of the near-degenerate ground states. Note that stim- 
ulated Raman adiabatic passage (STIRAP) is commonly used 
in atomic and molecular physics for the coherent transfer of 
quantum states 1 10], a task that is somewhat less general than 
the realization of arbitrary single-qubit operations. For an al- 
ternative proposal for single-qubit gates based on STIRAP see 
Ref. (H. 

In a real experiment, the adiabatic approximation will never 
be perfectly valid. Furthermore, spontaneous emission is an 
unavoidable error mechanism in any Raman system. In quan- 
tum dots, the interaction of the exciton with phonons can also 
be important, however it is possible to fabricate dots where 
spontaneous emission dominates all other sources of decoher- 



ence lT3l . In Ref. 0] the authors briefly discussed the condi- 
tions for adiabaticity and the effect of spontaneous emission 
from the excited state. In the present work we perform a more 
detailed study of these fundamental sources of error for the 
proposed gate protocol. In particular, we obtain quantitative 
results for the errors due to non-adiabaticity, and a simple for- 
mula for the errors due to spontaneous emission, namely that 
they are approximately equal to A7/A, where A is the rota- 
tion angle for the single-qubit rotation, 7 is the spontaneous 
emission rate and A is the detuning of the two Raman lasers 
from the excited state. We give evidence for this result both 
with a simple formal argument and by numerical computation. 

This paper is organized as follows. In section UU we de- 
scribe the protocol for realizing arbitrary single-qubit gates 
via Raman transitions proposed in Ref. 0]. In sectionlllllwe 
study the errors due to non-perfect adiabaticity (in the absence 
of spontaneous emission). In section HVl we study the errors 
caused by spontaneous emission. In section [V] we give our 
conclusions. 



II. GATE PROTOCOL 

In this section we describe the gate protocol proposed in 
Ref. 0. Consider a three level system composed of the two 
logical states of the qubit, |0) and |1), and an auxiliary excited 
state \X). The protocol relies on adiabatic Raman transitions 
in this A system, cf. Fig. [2 using two phase-locked laser 
pulses. The two laser frequencies are chosen such that they 
have the same detuning A, i.e. they satisfy the Raman reso- 
nance condition. 

The Hamiltonian in the interaction picture is: 



H = 




Oi(t)e ia 
fi a (t) 

n 2 (*) a 



(l) 



where a is the relative phase between the two real Rabi fre- 
quencies fti(t), Cl2(t). This Hamiltonian can be diagonalized 
straightforwardly. One introduces the following parameters: 
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FIG. 1: Scheme for single qubit gates analyzed in this paper. The 
low lying states |0) and |1) serve as qubit states. They are coupled to 
an excited state \X) via two laser beams with time-dependent Rabi 
frequencies Qi and fhe 1 ". The laser beams have the same constant 
detuning A from the excited state. 



The first eigenstate is time independent and com- 

pletely decoupled from the other two eigenstates. It has no 
contribution from the excited state \X). The second eigen- 
state \&2(t)) possesses only a small component of the excited 
state \X), as long as tl/A is small. The last state |$3(t)) is 
mainly composed of the excited state \X). 

An arbitrary unitary transformation can be realized adia- 
batically in the following way. Before the lasers are turned 
on (i.e. for (j> = 0), the initial qubit state can be expressed 
as a linear combination of the first two eigenstates. By apply- 
ing the two lasers, the Hamiltonian is then changed continu- 
ously. The adiabatic theorem states that, if the change of the 
Hamiltonian is sufficiently slow, the population in each instan- 
taneous eigenstate remains constant, only the relative phases 
of the eigenstates change. In the subspace formed by the first 
two eigenstates, one obtains the following transformation: 

[a,b]^ [a,be~^] (9) 

with 



z(t) = 4/n*(t) + 



i 



= |arctan^2^^ 



0(t) = arctan 



(3) 



(4) 



(5) 



The angle [3 is maintained constant by choosing the same 
envelope shape for the two pulses. One obtains the eigenval- 
ues : 



Ai(t) = with eigenvector 



l*i(*)> 



-e ta sin(/3) 
cos(/3) 




(6) 



M{t) = -2 Z(t) sin 2 (^(i)) with eigenvector 

— e ia cos(p) cos(^>) 
-sin(/3)cos(<£) 
sm(4>) 

X 3 (t) = 2 Z(t) cos 2 (0(i)) with eigenvector 

e la cos(/3) sin(</>) 
sin(/3) sin(0) 
cos(</>) 



(7) 



(8) 



A 2 



if 



\2(u)du, 



(10) 



where ti and t / denote the initial and final times respectively. 
From the point of view of the qubit basis spanned by the states 
|0) and |1), the resulting transformation has the form 



U 



-i/2A 2 3-n 



(ID 



where a is the vector of Pauli matrices, corresponding to a 
rotation through an angle A2 about an axis described by a unit 
vector n with components 



n x = cos(a) sin(2 (3) 
n y — — sin(a) sin(2 (3) 
n z = cos(2 (3). 



(12) 



We can now begin our detailed study of the corrections to this 
idealized description under realistic experimental conditions. 
We will start with errors due to non-perfect adiabaticity. 



HI. ERRORS DUE TO NON-ADIABATICITY 

A. Exact equations of motion 

The exact wave function can be expanded in terms of the 
previously defined instantaneous eigenstates, but with in gen- 
eral time-dependent coefficients: 



(u)du 



(13) 



Writing out the Schrodinger equation for this wave function, 
one obtains the following evolution equations for the coeffi- 



3 



cients: 



a k exp 



i / \ k {u)du 



(14) 



Y, a ^) ( $ fc(*) II CX P A »H rf 



Substituting the values of the scalar products (<bi(t)\\$j(t) 
according to the definition of the eigenstates, one finally has: 



di(t) 
aa(t) 
oa(t) 



= 



with 



-0(t) a 3 (t)P 23 (t) 
m a 2 (t)P 32 (t) 



(15) 



xmax 
7 



3.16 10 31.6 100 



FIG. 2: Determination of a; maa; , the maximum value of the ratio 
f2/A, as a function of the dimensionless parameter \ = At for a 7r 
rotation. The values for other rotation angles can be obtained from 
the fact that x max depends only on the ratio x/A. 



P32® = P23W = exp (i J (X 3 (v) ~ X 2 (v)) dv^j 

= exp (^2i jf Z(v)cfo^ . (16) 

The resolution of this system of differential equations allows 
to determine the error due to non-adiabaticity. 

For a given desired transformation U, we will define the er- 
ror as the maximum departure (in terms of overlap) of the real 
final state p(tf) from the ideal final state \4>ideai) — U\ip(ti)), 
where the maximization is over all initial states: 

E(U) = max [1 - (i>ideal\ P(tf) \lpideal)] ■ (17) 
V>(*i) 

The error can be expressed in terms of the complex coeffi- 
cients cti(tf). Since the coefficient cti always remains con- 
stant, a 3 (U) = and the differential equations are linear and 
homogeneous, it is in fact sufficient to solve the system of 
equations for one initial value of a 2 , say a 2 (ti) = 1, The re- 
sults for all possible initial states can then simply be obtained 
by multiplying the solution with the corresponding value of 
a 2 {U). 



To gain a better comprehension of the adiabatic approxima- 
tion, we make the following substitutions : for each function 
g of the time t, we write g(u) = g(ru), with the correspon- 
dence t = tu. Thus the new system evolves between the 
unitless time Ui = ti/r and Uf = tf/r. 

Introducing the dimensionless quantity 

X = At, (19) 

i.e. the product of the detuning and the gate time, we obtain 
the functions 

P 32 (u)=exp(-i X J u ^l + 4xl ax p(v)dv^j (20) 

and 

7/ \ %maxf{ u ) /TIN 

W = — jr—, (21) 

which appear in the dimensionless evolution equations: 

a 2 (u) = 4>(u) a 3 (u) P 32 (u) (22) 



B. Calculation of the error due to non-adiabaticity 

We now proceed to calculate the error due to non-perfect 
adiabaticity. To simplify the discussion, we will only con- 
sider laser pulse shapes /(£) that are approximate Gaussians 
of halfwidth r centered at t = 0, slightly modified such that 
the Rabi frequency is exactly zero at the initial and final times 

and tf). We introduce the ratio 

x(t) = ^-=x max f{t), (18) 

where x max is the maximal value of the ratio CI/ A, i.e. we 
have normalized / such that /(0) = 1, 



a 3 (u) = -4>(u) a 2 (u) P 32 {u) (23) 

We thus obtain a new system of differential equations de- 
pending on the two dimensionless parameters \ et x max- We 
are interested in the dependence of the solutions on the two 
parameters. First it appears that the greater \ is, the faster the 
term P 32 is oscillating, and thus the less the population |a,3 1 2 
is important. We can also obtain a reduction of 1 03 1 by reduc- 
ing (p. This can be done by decreasing x max . This preliminary 
analysis suggests that the error decreases with \ and increases 
with x max . 

In the following, we will study the gate error as a function 
of the rotation angle A and the dimensionless parameter x- 
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FIG. 3: Gate error as a function of \ for different rotation angles. 



The quantity x m ax is then not an independent variable, but 
is determined by these two parameters in the following way. 
After simplifications and substitutions, Eq. (II Oi becomes : 



2A = X J f V 1 + 4 X ™aJ 2 ( U ) ~ljdu = X 9{Xmax). 

(24) 

Recall that f(u) is essentially a Gaussian with halfwidth 1 
(apart from a small modification at the boundaries of the time 
interval), normalized such that its maximum value is equal to 
one. This equation gives x max as an implicit function of the 
ratio A/x- As g is an increasing function, and is a bijection 
from [0, +oo] to [0, +oo], we obtain a one to one correspon- 
dence between x max and \ for a given A. Fig. |2]shows that 
Xmax is a decreasing function of x- 

For a given rotation angle A, we finally obtain a system 
depending only on the parameter x- By solving it numerically, 
we obtain an estimation of the adiabatic error as a function 
of x an d of A. Note that the error does not depend on the 
axis of rotation. Indeed, the choice of the axis of rotation 
only determines the relation of the basis of logical states to the 
basis of adiabatic eigenstates, and since the error is obtained 
by maximizing over all initial states, this choice has no effect 
on its value. The error as a function of x for different values 
of the angle A is shown in Fig. [3] For each value of A, we 
observe two characteristic regimes. The first one, where x is 
small, is characterized by a damped oscillatory behaviour as a 
function of x, leading also to a non-monotonous variation of 
the error with the rotation angle. In the second one, the error 
decreases continuously with x, and greater values of A lead to 
larger errors. 

In order to minimize the error due to non-adiabaticity, 
X = At should thus be as large as possible. However, it 
should be kept in mind that r corresponds to the time of the 
gate and thus has to be much shorter than the decoherence 
time of the qubit states, in order to limit errors due to decoher- 
ence. In principle one can be in the adiabatic regime even for 
very short gate times, provided that the detuning A is made 
large enough. However, this requires an increase in the laser 
amplitude CI, in order to still achieve the same rotation an- 
gle A. This relationship is contained in Eq. i24\ and Fig. |2j 
which show that for fixed x the rotation angle is determined 



by x max , i.e. the ratio of CI and A at maximum laser intensity. 
If one decreases r and increases A, keeping x an d thus the 
level of error constant, one therefore has to increase CI by the 
same factor as A, in order to keep the rotation angle constant. 
Since the laser intensity cannot be made arbitrarily large, this 
imposes an upper bound on A, and thus a lower bound on r. 
From Fig. |3]one can see that for x > 15, the error is signifi- 
cantly less than 10~ 4 , which should be small enough for fault 
tolerant quantum computation QJ . 

IV. ERRORS DUE TO SPONTANEOUS EMISSION 
A. Estimate of error based on population transfer 

In this section we will investigate the error introduced to the 
Raman single-qubit gates by the finite lifetime of the excited 
state | A) due to spontaneous emission. We will begin with a 
fairly simple argument that gives the correct behaviour for the 
error, before presenting more precise numerical calculations 
in the next subsection. 

As described before, the two adiabatic basis states that are 
significantly populated during the gate operation are 
and |$2(i))- The state has no contribution from the 

excited state and is thus unaffected by spontaneous emission. 
On the other hand, the state |$2(0) has a component in the 
excited state |A). From Eq. (Q, we can see that the popu- 
lation in the unstable state \X) is therefore a\(t) sin 2 (0(i)). 
Defining 7 as the spontaneous decay rate of the state |A), we 
can then estimate the population <5 transferred by spontaneous 
emission during the gate operation as follows: 

6 = S 7 aj(t) sin 2 (<p{t))dt (25) 

Let us assume that we are in the adiabatic regime (x > 15), 
and that the overall error due to spontaneous emission is small. 
The first of these conditions implies that x max and thus <j> is 
small, cf. Fig. [2] The second one implies that 02 is nearly 
constant. Eq. (I25> can then be simplified to 

S = lTalxl iax f 2 (u)du, (26) 

J Ui 

where we have again introduced the dimensionless function 
f(u) defined above. Furthermore, in the same regime, Eq. 
(I24> can be simplified to 

r u f _ A 
a&a* / A«)d« = -■ (27) 

Choosing a| = 1 in order to obtain an upper bound, we find 
the following expression for the total transferred population 
due to spontaneous emission: 

5 = Al (28) 

The error induced by spontaneous emission is twofold; on the 
one hand, a new distribution of the populations | of | , and on 
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the other hand, a dephasing between the two qubit basis states. 
The transferred population S provides an estimate for the gate 
error due to spontaneous emission. It may seem surprising 
that 6 does not depend on the gate time r, even though for 
longer gate times the component of the system in the excited 
state has more time to decay. The reason for this is that for 
the same rotation angle shorter gate times require larger pop- 
ulations in the excited state, and the two effects cancel out ex- 
actly, at least within the framework of the above estimate. We 
are now going to use numerical computation to obtain more 
precise results on the gate errors. 



Population 
1 < 



0.8 
0.6 
0.4 
0.2 



-0.02 




(a) 



0.02 0.04 



t (ns) 



B. Master equation 

In order to study the effects of spontaneous emission on the 
Raman gate protocol in detail, we use the master equation for- 
malism. In the present section, we assume for simplicity that 
the spontaneous emission can only occur from the state \X) 
toward the two qubit states |0) and |1), and not to any addi- 
tional states. To take this decay into account, we introduce the 
Lindblad operators : 



Population 
1 1 



0.8 
0.6 
0.4 
0.2 



-0.02 




(b) 



0.02 0.04 



t (ns) 



L, 



(29) 



with <j i = \i) (X\ (for i — 0, 1). The constants ji are the 
decay rates from \X) towards the states \i). The total decay 
rate is thus given by 7 = 71 +72. 

The master equation {jl is: 



dp 



-i[H,p]+J2(2L iP Lf -{L+L l:P }) 



(30) 



For quantum gate operations, the initial state \ip) is always a 
linear combination of the two logical states a |0) + b |1) with 
\a\ 2 + \b\ 2 = l. 

The master equation corresponds to a set of coupled dif- 
ferential equations for the elements of the density matrix that 
can be solved numerically. This allows us to determine the 
gate error defined in Eq. (I17> . 



C. Example - Populations and Purity 

In this subsection, we describe the numerical results for a 
particular case in detail, comparing the situations with and 
without spontaneous emission. This is intended to serve as 
an introduction and illustration for our more general results 
presented in the next subsection. 

We consider a rotation by n along the x axis, correspond- 
ing to a — and j3 = 7r/4. . The initial state of the sys- 
tem is the state |0). The ideal final state is the state |1). We 
choose the values r = 0.01 ns, A = 1 meV« 1500 ns -1 and 
7i = 72 = 20 ns^ 1 . While these values could apply to con- 
ceivable experiments with quantum dots fl3ll . we have also 
chosen them such that the relevant effects are clearly visible. 
For these values \ = 15. The adiabatic approximation is thus 
very well satisfied, cf. Fig. [3] 



FIG. 4: Evolution of the populations during a n rotation of the state 
1 0) around the x axis into the state 1 1 ) : (a) without spontaneous emis- 
sion; (b) with spontaneous emission. 



Fig. |4ja) shows the time evolution of the populations in the 
states |0), |1) and \X) without spontaneous emission. As ex- 
pected, the populations in the states |0) and |1) are exchanged. 
Moreover, the population in \X) is zero at the end of the oper- 
ation. One can describe the evolution of the system in detail as 
follows. The initial state |0) can be expressed in the adiabatic 
basis as: 



10) 



1 



V2 



(l*i(*0> + l*a(<0» 



The final state |1) can be written : 



|1) 



V2 



(31) 



(32) 



with |$i(t/)) - and |$ 2 (t/)) = |* a (fi)>. 

The transformation is adiabatic : expressed in the basis 
(|$i(£)),|$ 2 (*)))> the state of the system is : 



w» = -> 



iA(i) 



(33) 



While the lasers are on, the phase A(i) grows and thus the 
state |0) is continuously transformed into the state |1). The 
state I ) ) has a component of order x max f(t) in the state 
\X), the population in the excited state therefore grows until 
the maximum of the light intensity is reached, and then returns 
to zero. Fig. |4jb) shows the same populations in the presence 
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FIG. 5: Evolution of the purity of the state, Q(t), without and with 
spontaneous emission. 
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FIG. 7: Gate error in the presence of spontaneous emission as 
a function of the gate time r. The spontaneous decay rates are 
7i = 72 = 5 ns -1 , the detuning A = 1 meV. The variation of 
r between and 20 ps corresponds to the dimensionless parameter 
X = At varying from to 30. 



the presence of spontaneous emission causes transitions be- 
tween the states |$2(i)) and |$3(f)), which populate the lat- 
ter. This population decreases on the timescale of the radia- 
tive lifetime, since | <3>3 (i)) contains predominantly the excited 
state I JO. 



D. General results on errors due to spontaneous emission 



FIG. 6: Time evolution of the population in the adiabatic basis state 
|<3?3(t)) (whose dominant component is in the excited state \X)) 
without and with spontaneous emission. 



of spontaneous emission. It appears clearly that the rotation 
is no longer perfect in this case. The populations in the initial 
state |0) and in the excited state \X) are no longer zero at the 
moment when the light is turned off. Of course the remaining 
population in \X) will decay towards the states |0) and |1) on 
the larger timescale set by ji 2 . 

The departure from the perfect rotation can be further vi- 
sualized by analyzing the purity of the system density matrix 
pit), i.e. by studying the quantity 0(i) = Trp(i) 2 . Fig. [5] 
shows the evolution of with and without spontaneous emis- 
sion. As expected, remains constant in the absence of spon- 
taneous emission. In the presence of spontaneous emission, 
the purity decreases considerably. This decrease is particu- 
larly strong around t = 0, i.e. around the maximum of the 
laser intensity. It is at this point that the population in the state 
\X) becomes the largest. The spontaneous emission from \X) 
causes the state to become mixed. 

Fig. |6] shows the time evolution of the population in the 
state | $3(0) of Eq. (|8}. One sees that in the absence of spon- 
taneous emission the adiabatic approximation is well justified, 
the population in the state remains very small. Its departure 
from zero corresponds to the error due to non-perfect adia- 
baticity discussed in the previous section. On the other hand, 



In this subsection we will present more general results on 
the error due to spontaneous emission. In particular we want 
to test the estimate made in subsection IIV Al Fig. shows 
the error for a rotation by it as a function of the gate time r in 
the presence of a spontaneous emission. This graph should be 
compared to Fig. [3J which shows the same quantity in the ab- 
sence of spontaneous emission. For short times the behaviour 
is very similar, showing the same damped oscillatory charac- 
ter. In this regime, the error is dominated by non-perfect adi- 
abaticity. For longer gate times, there is a clear difference. In 
the presence of spontaneous emission, the error does not fall 
below a certain minimal value and is virtually independent of 
t. This is in good correspondence with the prediction made in 
subsection llV Al 

In order to test the above estimate more systematically, we 
have performed calculations varying 7 and A for a fixed value 
of the gate time r = 13.3 ps. We restrict ourselves to values of 
A > 1 meV, corresponding to \ > 20, in order to make sure 
that the error due to non-adiabaticity is negligible, cf. Fig. [3] 
Furthermore we focus on the regime where the overall gate 
error is at most at the percent level, since this is the relevant 
regime for quantum computing. The total spontaneous decay 
rate is 7 = 71+72- For simplicity we have again chosen 
7i = 72- 

Fig. [8]shows the behaviour of the gate error of a tt rotation 
as a function of the spontaneous decay rate 7 for four different 
values of the detuning A. One sees that the error is linear in 
7 with a high degree of accuracy. Similarly, Fig. [9] shows 
the gate error as a function of A for three different values of 
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FIG. 8: Gate error for a n rotation as a function of the spontaneous 
decay rate 7 for fixed values of the detuning A = 1, 2, 4, 8 meV 
respectively (from top to bottom graph). 




2 3 4 5 6 



A (meV) 



FIG. 9: Gate error of a n rotation as a function of detuning A for 
fixed values of the spontaneous decay rate 7 = 2, 4, 6 ns _1 respec- 
tively (from bottom to top graph). The curves are fits to a 1/ A be- 
haviour. 



7. One sees that the results are fitted extremely well by a 1/ A 
behaviour. The proportionality of the error to 7/ A is thus seen 
to be very well obeyed. To assess the accuracy of the above 
estimate concerning the absolute size of the error, Fig. ^| 
compares the errors obtained numerically for a it rotation to 
the estimated error of irj/A. One sees that the approximation 
works very well in the considered regime. As expected, it 
tends to work somewhat less well for increasing values of 7 
and decreasing values of A, i.e. increasing overall size of the 
error, cf. subsection llV Al 



V. CONCLUSIONS 

The results obtained in the present paper give quantitative 
information for the implementation of quantum computing us- 
ing the considered gate protocol. The results of sectionliTflon 
the errors due to non-perfect adiabaticity make it possible to 
determine the maximum allowable gate speed for any desired 
level of error. Our analysis also shows that it is in principle 
possible to perform gates in the adiabatic regime even for very 



ratio 




FIG. 10: Ratio of the (numerically obtained) exact value of the er- 
ror for a 7T rotation to the value of ^7/ A given by the estimate of 
subsection llV Al for a range of values of 7 and A. 



short gate times by increasing the detuning A. However, this 
requires a corresponding increase in laser power. 

The results of section ffvl quantify the errors due to spon- 
taneous emission, whose presence is unavoidable in any gate 
scheme based on Raman transitions. In the context of quan- 
tum computing with spins in quantum dots, the present analy- 
sis complements the results of Refs. 1161 [1411 on quantum gate 
errors due to phonon-induced dephasing. Phonon-related er- 
rors can be made essentially arbitrarily small by making the 
gate operation slower. The basic reason for this is that the 
speed of the operation determines the energy that is available 
for the creation of phonons (since the system does not decay 
from the excited state during the dephasing). The slower the 
operation, the less energy is available, restricting the available 
state space for phonon creation. Unfortunately there is no cor- 
responding energy constraint for spontaneous emission, since 
the energy for photon creation is provided by the decay of the 
emitter to one of the low-lying states. Our results show ex- 
plicitly that slowing down the gate operation is not helpful to 
reduce errors in the present context, see Fig. [7] 

On the other hand, our analysis also shows that it is possi- 
ble to choose long gate times even in the presence of sponta- 
neous emission, without significantly changing the size of the 
error due to the decay. Relatively long gate times can be ad- 
vantageous because they allow greater frequency selectivity in 
schemes based on spectral addressing of different qubits. Of 
course, the gate time always has to be much shorter than the 
decoherence time of superpositions of the qubit states. 
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